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ABSTRACT 


Flood  frequency  analysis  for  gauged  and  ungauged  sites  in  Barak  river  basin  in  south  Assam  are  conducted  using 
peak  discharge  and  peak  gauge  level  for  the  sites.  The  best  fit  distribution  for  the  regions  were  selected  based  on  the  results 
of  ZDIST  statistics  and  L-moments  ratio  diagram  from  the  five  extreme  value  distributions  GPD,  GLO,  GEV,  P3  and  LP3. 
Regional  flood  frequency  relationships  were  established  using  the  selected  distribution  along  with  development  of  growth 
factors  for  the  sites.  A  regional  flood  relationship  with  the  physical  characteristics  for  the  sites  was  established  using 
multiple  regressions  for  the  estimation  of  flood  at  the  ungauged  sites  with  the  validity  testing  of  the  developed  equation. 

KEYWORDS:  Growth  Factors,  L-Moments,  L-Moments  Ratio  Diagram,  ZDIST  Statistics 


The  observed  flood  flow  in  terms  of  discharges  and  flood  levels  for  rivers  and  streams  in  a  drainage  area  can  be 
used  in  various  engineering  works  like  design  of  dams  and  spillways,  flood  control  works,  design  of  bridges  and 
construction  of  roads  and  railways  and  forecasting  of  flood  warning  system.  Flood  frequency  analysis  is  to  determine  how 
large  a  flood  of  a  certain  magnitude  in  terms  of  discharge  or  flood  level  can  occur  at  a  given  site  during  the  desired  return 
periods.  This  estimated  flood  magnitude  can  give  the  most  vital  information  during  the  planning  of  many  engineering  and 
flood  related  works.  Frequency  analysis  at  a  single  site  can  give  this  desired  information  if  the  observed  data  from  the  site 
are  of  sufficient  length  in  comparison  to  the  desired  return  periods.  In  case  of  the  sites  where  the  available  data  are  of  very 
short  period  or  data  from  the  sites  are  not  available,  at  site  method  of  frequency  analysis  can  not  be  able  to  give  valid  and 
reliable  quantile  estimation.  Regional  based  frequency  analysis  is  the  alternative  method  of  frequency  analysis  where  the 
available  data  from  the  different  sites  either  long  or  short  periods  are  pooled  together  so  as  to  augment  the  site's  data  for 
the  estimation  to  a  site  in  a  region.  L-moments  based  regional  flood  frequency  analysis  approach  which  is  incorporated 
with  the  identification  of  discordant  sites  along  with  the  testing  of  regional  homogeneity  is  now  the  most  commonly  used 
method  everywhere.  L-moments  are  the  linear  combination  of  probability  weighted  moments  (Hosking  1986,  1990). 
They  are  less  biased  for  the  estimation  of  parameters  of  frequency  distributions  and  less  effective  in  the  process  of 
frequency  analysis  during  the  presence  of  outliers  in  the  data. 

L-moment  approach  has  been  used  for  regional  flood  frequency  analysis  (Bardrelin  GH.  Hassan  et  al.  2011  in 
China,  Hussein  Maleki-Nezhad  et  al.  2011  in  Iran,  Sattar  Chavoshi  Borujeni  et  al.  2009  in  Malaysia,  Ramin  Rostami  2013 
in  Iran,  Yeo  H.  Lim  et  al.  2009  in  USA  ,  Anil  Kumar  et  al.2012  and  Rakesh  Kumar  et  al.,  2005  in  India).  In  this  study 
L-moments  based  regional  flood  frequency  analysis  for  the  gauged  and  ungauged  sites  for  Barak  Valley  of  Assam,  India 
has  been  conducted  using  both  annual  peak  discharge  and  annual  gauge  height.  So  far  no  such  study  has  been  found  for 
flood  frequency  analysis  for  the  ungauged  sites  in  this  valley.  Hence,  this  study  has  been  carried  out  for  the  estimation  of 
T-years  flood  in  terms  of  peak  discharge  and  gauge  height  for  both  gauged  and  ungauged  sites  in  this  valley. 


INTRODUCTION 
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The  objective  of  this  study  is  to  apply  L-moment  based  approach  in  regional  extreme  flood  frequency  analysis  for 
the  gauged  and  ungauged  sites  of  Barak  valley  of  India  and  to  establish  regional  mean  flood  relationship  with  physical 
characteristics  of  the  sites  using  multiple  regressions  for  the  estimation  of  flood  at  the  ungauged  sites. 

STUDY  AREA 

This  study  area  lies  within  24°  N  to  26°  N  and  92°  E  and  94°E  in  the  southern  part  of  Assam,  North  East  India. 
This  valley  comprises  of  three  districts  Cachar,  Hailakandi  and  Karimganj  of  Assam.  The  annual  rainfall  of  this  valley  is 
2500-3000  mm.  Most  of  the  floodplains  in  Barak  valley  are  in  the  low-lying  areas  that  get  inundated  during  the  monsoon 
months  of  June-September  by  the  flood  waters  of  Barak  and  its  tributaries.  The  Barak  river  originates  from  the  hills  of 
Manipur  and  flows  touching  Mizoram  State  before  it  enters  the  plains  of  Assam  at  Lakhipur.  After  flowing  through 
Silchar  it  bifurcates  near  Badarpur  into  Surma  river  and  Kushiyara  river  and  then  it  enters  into  Bangladesh.  The  principal 
tributaries  of  the  Barak  in  India  are  Jiri,  Dhaleshwari  (Tlawng),  Singla,  Longai,  Madhura,  Sonai  (Tuirial),  Rukni,  Jatinga 
and  Katakhal.  The  annual  peak  discharge  data  from  7  gauged  sites  and  peak  gauge  level  from  8  sites  for  the  period  from 
1996  to  2010  spread  over  the  valley  are  collected  from  Central  Water  Commission  of  India,  Shillong.  The  map  of  the 
study  area  is  given  in  Figure  1  and  the  summary  of  the  gauging  sites  are  listed  in  Table  1 .  The  steps  for  L-moment  based 
frequency  analysis  approach  with  the  index  flood  procedure  are  given  briefly  in  the  proceeding  sections. 
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Figure  1:  Map  of  the  Study  Area  (Enclosed  in  Red  Line) 
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Table  1:  Details  of  the  Gauged  Sites  with  the  Availability  of  the  Data 


Station 

River 

Latitude 

Longitude 

Available  data 

Catchment 
Area  (km2) 

Discharge 

Gauge 

l.A.P.Ghat 

Barak 

24°49'58"N 

92°47'30"  E 

do 

do 

18721 

2.Fulertal 

Barak 

24°47'19"N 

93°01'08"E 

do 

do 

14450 

3.Gumra 

Gumra 

25°00'41"N 

92°30'  35"E 

do 

do 

2800 

4.Dholai 

Rukni 

24°35'10"N 

92°50'32"  E 

do 

do 

562 

5.B.P.Ghat 

Barak 

24°52'32"N 

92°35'00"  E 

do 

do 

24216 

6.Karimganj 

Kushiyara 

24°52'N 

92°21'  E 

do 

7.Fakirabazar 

Longai 

24°51'06"N 

92°20'43"  E 

do 

do 

1108 

8.Matizuri 

Katakhal 

24°38'53"N 

92°36'  29"E 

do 

do 

7770 

,  C/-1JC  j-aK;-a3..i.i.i.rj-L] 
-lJO-eJCw-s)  (W-O 


L-MOMENT  APPROACH 

L-moments  are  linear  combinations  of  probability  weighted  moments  (PWM).  The  probability  weighted  moments 
are  calculated  from  the  ranked  observations  (Xj).  Greenwood  et  al.  (1979)  defined  the  probability  weighted  moments  as 

The  first  four  L-moments  in  terms  of  probability  weighted  moments  are  written  by 

*l  =  ft 
2S  =   2ft  -  ft 

2S=  6ft  -  6ft  +ft 

A4=  20ft -30  ft  +  12ft -ft  (Id) 


(1) 

(la) 

(lb) 
(lc) 


where  Xt  is  the  L-mean  which  measures  the  central  tendency,  Jl :  is  the  L-standard  deviation  and  measures  the 
dispersion.  The  dimensionless  L-moment  ratios  defined  by  Hosking  (1990)  are  written  as 


L—  vaaf f  cis-nt  af  va-ritmcBiT  =  -r 


I  -  Skewness:  zt  =  — 


L  —  Kvrtasisnj,  =  — 


(le) 
(If) 
dg) 


The  regional  weighted  values  of  A\  X2,  t-,  t3  and  t4  for  N  sites  in  a  region  can  be  obtained  from  the  following 


relations 


c  S.-=1Cniii]  E.-^tfiiiji]  S,=1tnL^i3  _  %i=f*i*3£  _  S,-=i"iT*i- 
R  _   iZt  .  r  R  Lit  .  „s  Lit  _.  r  Lit  s  _   Lit  

Sj=/«l3  2J=1^il  2j=1W  2>ltni3 


Where  i  =  {  site  and  N  =  number  of  sites.  These  regional  values  of  weighted  L-moments  and  their  ratios  are  the 
basic  parameters  for  L-moments  based  regional  frequency  analysis. 
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DISCORDANCY  MEASURE 

The  aim  of  the  regional  frequency  analysis  is  to  select  the  best  fit  frequency  distribution  applicable  to  all  the  sites 
over  the  region.  So,  the  collected  data  from  the  sites  should  be  from  the  same  frequency  distribution.  Discordancy  measure 
is  to  screen  out  data  from  the  unusual  sites  whose  sample  L-moments  are  markedly  different  from  other  sites.  It  is 
measured  in  terms  of  L-moments  (Hosking  &  Wallis,  1993)  and  is  defined  as 

A=  mFA_l(m£-  u  )  (2) 

Where  u£  =  (t,  )  T  is  a  vector  containing,  i,  t3  and  t,  values  of  a  site  i  ,  the  superscript  3"  denotes 
transpose  of  a  matrix  or  vector;  it  =  ^  £k£  represents  the  group  average  and  A~L  is  the  inverse  of  the  covariance 
matrix  A  of  u^.The  elements  in  A)'1-  are  given  by  the  following  relation 

A=  I  J&iCui-  u)CC«i-  (3) 

Where  N  is  the  number  of  sites  in  the  region.  The  station  ;  is  declared  as  discordant,  if  its  (  D £  ]  value  obtained 
using  equation  (2)  is  greater  than  the  critical  value.  The  critical  values  of  (  D^}  depend  upon  the  number  of  sites  in  the 
region.  For  a  region  having  a  number  of  sites  >15,  the  critical  value  is  taken  as  3  however,  if  the  number  of  sites  is  <  15  the 
critical  value  can  be  obtained  from  the  table  given  by  Hosking  and  Wallis  (1993). 

The  L-moments  and  their  ratios  with  the  discordancy  (D^J  for  the  7  sites  using  annual  peak  discharge  data  and 
8  sites  using  annual  peak  gauge  level  data  are  given  in  Table  2  and  Table  3.  It  is  observed  that  there  is  no  discordance  site 
available  for  both  the  data.  The  regional  L-moments  and  their  ratios  for  these  sites  are  in  Table  4. 


Table  2:  L-Moments  and  Discordancy  Values  for  the  Sites  (Discharge) 


Station 

n 

T 

T4 

^critical 

D, 

l.A.P.Ghat 

13 

20.60915 

0.487051 

0.023633 

-0.11091 

0.191232 

2.140forN  =  8 

0.6189 

2.Fulertal 

13 

24.83423 

0.50981 

0.020529 

0.108812 

0.042471 

0.7323 

3.Gumra 

13 

17.65308 

0.17218 

0.009754 

-0.0429 

0.031006 

1.4521 

4.Dholai 

13 

27.23077 

0.328203 

0.012053 

-0.33857 

0.376962 

1.7765 

5.B.P.Ghat 

13 

17.51654 

0.322052 

0.018386 

-0.0317 

0.276041 

0.7103 

6.Karimganj 

15 

15.87467 

0.245905 

0.01549 

0.03803 

0.175753 

0.2793 

7.Fakirabazar 

13 

17.22039 

0.256218 

0.014879 

0.050554 

0.056539 

0.4048 

8.Matizuri 

13 

21.81808 

0.548269 

0.025129 

-0.41457 

0.08661 

2.0254 

Table  3:  L-Moments  and  Discordancy  Values  for  the  Sites  (Gauge  Level) 


Station 

n 

To 

T4 

^critical 

A 

l.A.P.Ghat 

13 

20.60915 

0.487051 

0.023633 

-0.11091 

0.191232 

2.140forN  =  8 

0.6189 

2.Fulertal 

13 

24.83423 

0.50981 

0.020529 

0.108812 

0.042471 

0.7323 

3.Gumra 

13 

17.65308 

0.17218 

0.009754 

-0.0429 

0.031006 

1.4521 

4.Dholai 

13 

27.23077 

0.328203 

0.012053 

-0.33857 

0.376962 

1.7765 

5.B.P.Ghat 

13 

17.51654 

0.322052 

0.018386 

-0.0317 

0.276041 

0.7103 

6.Karimganj 

15 

15.87467 

0.245905 

0.01549 

0.03803 

0.175753 

0.2793 

Impact  Factor  (JCC):  5.7179 
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Table  3:  contd., 


7.Fakirabazar 

13 

17.22039 

0.256218 

0.014879 

0.050554 

0.056539 

0.4048 

8.Matizuri 

13 

21.81808 

0.548269 

0.025129 

-0.41457 

0.08661 

2.0254 

Table  4  :  Regional  L-Moments  and  their  Ratios 


Regional  L-Moments 
and  Their  Ratios 

Discharge 

Gauge 

3% 

2114.0166 

20.260273 

-  R 

239.34528 

0.356582 

0.1324421 

0.017443 

0.0581456 

-.0.090192 

■  T 

0.1521240 

0.1549763 

HETEROGENEITY  MEASURE  (HOSKING  and  WALLIS,  1997) 

The  heterogeneity  measure  estimates  the  degree  of  heterogeneity  in  a  group  of  sites  and  determines  whether  the 
group  of  sites  can  be  reasonably  treated  as  a  homogeneous  region.  It  compares  the  inter-site  variations  in  sample 
L-moments  for  the  group  of  sites  with  what  would  be  expected  of  a  homogeneous  region.  A  statistic  based  on  the  weighted 
variance  of  the  L-coefficient  of  variation  (L-cv)  is  derived  and  heterogeneity  measure  is  defined  as 

H  =  ^  (4) 
where  Y  =  weighted  standard  deviation  of  L-cv  for  the  observed  data  and  is  given  by 

Here,  «j  =  record  length  for  the  site  ;  t  =  the  L-Cv  for  the  site     i  =  1,2  ,3....   =  the  weighted  group 

mean  L-Cv  ;  *V=  no.  of  sites  in  the  region  and  nL.and  DL,are  the  mean  and  standard  deviation  respectively  for  the  weighted 
standard  deviations  for  the  500  sets  of  simulated  homogeneous  regions  generated  by  Monte  Carlo  simulation  using  a 
4-parameter  Kappa  distribution  (Hosking,  1997). 

The  quantile  function  of  the  distribution  is  given  by 

The  parameters  of  this  distribution  f ,  a,  K  and  h  are  estimated  using  their  relationship  with  regional  A2R,  tr,  Tj" 
and  t4s  values;  F  is  the  probability  and  x(F)  is  the  flood  quantiles.  The  flood  quantiles  for  the  different  sites  are  simulated 
using  the  quantile  function,  given  in  equation  (6)  from  the  data. 

The  criteria  for  deciding  the  heterogeneity  of  a  region  is  as  follows 

•  if  H  <  1,  the  region  is  acceptably  homogeneous, 

•  if  1  <  H  <  2  ,  the  region  is  possibly  heterogeneous, 
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•     if  H  >  2  ,  the  region  is  definitely  heterogeneous. 

The  parameters  for  Kappa  distribution  f ,  a,  K  and  h  along  with  the  values  of  V,  m,,  av  and  H  for  the  7  and  8 
sites  are  given  in  Table  5.  The  H  values  for  the  7  and  8  sites  are  0.2515  and  -0.4088  respectively.  Hence  this  valley  as 
whole  is  homogeneous  for  both  the  data. 

Table  5:  Parameters  of  Kappa  Distribution  f,  a.  K  and  h  Along  with  the  Values  of  V,  juL„  fl,  and  H 


Parameters 

For  Annual  Peak  Dischare 

For  Annual  Peak  Gauge 

fi 

-0.564199 

-0.422399 

K 

0.157499 

0.004699 

a 

0.029254 

0.004415 

c 

0.2088678 

0.201074 

Homogeneity  Measure 

V 

0.03883 

0.005037 

0.03858 

0.005423 

0.00994 

0.000944 

H 

0.2515 

-0.4088 

GOODNESS  OF  FIT  MEASURE 
Statistics 

The  goodness  of  fit  criterion  for  each  of  the  candidate  distributions  is  carried  out  by  comparing  the  L-moments  of 
the  candidate  distributions  to  the  regional  average  L-moments  statistics  derived  from  the  regional  data.  For  each  of  the 
distributions,  the  goodness  of  fit  measure  (Hosking  and  Wallis,  1997)  is  given  by 

2D1ST  —      *         ~    <.   +  ^ 

where  t4s  =  weighted  L-kurtosis  obtained  from  the  data  of  a  given  region. 

t^;st  _  l_  kurtosis  of  the  fitted  distributions  which  are  derived  from  their  relationship  of  regional  L-kurtosis 
T4sirr/'  =  weighted  L-kurtosis  obtained  from  the  simulated  data. 

ffj  =  — '  ^   1    measures  bias  of  weighted  t4s  with  weighted  Tj^uj*  of  the  simulated  regions  and  tr4 

is  computed  using  the  relation 

i 

°4  =  kr^-r,  *  -        -  xsirK&  }f  (8) 

Here,  iV^j^the  number  of  simulated  regional  data  sets  generated  using  a  4-parameter  Kappa  distribution. 
If  l^'57!  obtained  by  using  equation  (12)  is  close  to  zero  the  fit  is  considered  to  be  adequate  and  reasonable  if 
|^s;st|  ^  1,54,  If  all  the  distributions  have  their  |Z£'i'ST|  <  1,64  ,  the  one  which  is  closest  to  zero  is  selected  as  the  best  fit 
distribution  for  a  region.  The  l^^53"!  values  for  both  the  data  are  given  in  Table  6.  It  is  observed  that  the  best  fit 
distribution  for  both  the  data  are  found  as  LN3. 
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Table  6:  Selection  of  Best  Fit  Distribution  Using  \Zmsi\  Statistics 


Distribution 

|ZDIST| 

Statistics 

L-Moment 

Ratio 
Diagram 

IZDISTI 
Statistics 

L-Moment 

Ratio 
Diagram 

GPD 

0.4623 

0.4948 

GLO 

0.1311 

0.1578 

GEV 

0.0752 

0.0931 

P3 

0.0853 

0.0649 

LN3 

0.0440 

LN3 

0.0098 

LN3 

L-MOMENTS  RATIO  DIAGRAM 

L-Moment  ratio  diagram  (Kroll  &  Vogel,  2002  and  Kumar  et  al,  2003)  is  another  tool  commonly  used  for  the 
selection  of  the  best  distribution  by  comparing  its  closeness  to  the  regional  L-skewness  and  L-kurtosis.  The  L-kurtosis  (  t4) 
for  each  of  the  distribution  are  derived  from  their  relationship  of  L-skewness  (t2}  for  the  sites.  The  L-kurtosis  (14)  derived 
for  each  of  the  distribution  are  plotted  against  L-skewness  (t,]  for  the  sites  and  the  distribution  which  is  closest  to  the 
regional  values  of  (tj'5',  T4?')  is  selected  as  the  best  fit  distribution. 

The  L- moments  ratio  diagram  for  both  the  data  are  given  in  Figure  2  and  Figure  3.  In  this  test  also  the  best  fit 
distribution  for  both  the  data  can  be  taken  as  LN3.. 

■  (8477  ■  0L0    ■  OEV    ■  P3  ■  LN3 — ■  T3-T4 — 

0.377   

0.277   


-0.022     1  1  1  1  1  1— 

-0.163  -0.063  0.036  0.136  0.236  0.336 

Tau3 

Figure  2:  L-Moment  Ratio  Diagram  (Discharge) 


www.tjprc.org 


editor@tjprc.org 


14 


Thokchom  Arti  Devi  &  Parthasarathi  Choudhury 


■   GPD    ■  GLO    ■  GEV    ■  P3       ■   LN3     ■  T3-T4 
0.389       |    1    1    1  


0.289 


T 

■3        0  189 

u 

4 


0.089 


-0.010 


-0.414 


0.185 


Figure  3:  L-Moment  Ratio  Diagram  (Gauge  Level) 

INDEX  FLOOD  PROCEDURE 

The  index  flood  procedure  assumes  that  in  a  homogeneous  region  the  frequency  distribution  at  the  sites  are 
approximately  identical  apart  from  a  scale  factor  (Dalrymple,  1960).  For  a  homogeneous  region  having  N  sites,  with  site  i 
having  sample  size  «,  and  observed  data        the  following  relationship  given  in  equation  (14)  holds  good  for  a  site  i  in  the 

homogeneous  region 

Qi{F}=  Vi?(F}  (9) 

where  j  =  1,2  ,3  rv     and    i=  1,2  ,3  N.;       00  ,  0  <  F  <  1,  is  the  quantile  function  of  the 

frequency  distribution  at  site  i.  \i;  is  a  site-dependent  scale  factor,  called  the  index  flood  and  is  often  taken  as  the  at  site 
mean;  tf{F)  =  regional  growth  factor  (Dalrymple  ,1960) 

If,  p[  be  the  scale  factor  or  at  site  mean  for  site  i,  then  the  dimensionless  rescaled  data  for  site  i  can  be  obtained  by  dividing 
the  observed  data      by  the  at  site  mean     =  1± 


qi!  ~    R  " 


(10) 


Using  Equation  (15)  a  new  set  of  rescaled  flood  data  for  each  of  the  stations  can  be  obtained.  The  regional 
L-moments  and  their  ratios        ,  A-jj   and  jw_a  of  the  rescaled  data  can  be  determined  as  done  in  case  of  the  observed 

data  or  it  can  be  taken  as  X  J  =  1 .0,  iR  =  ^  =  ^4  =        as  \u*±  1 .0  and         =  ^4  =  ^4  =  T, * 

The  parameters  of  the  selected  distribution  can  be  estimated  from  their  relationships  of  L-moments  and  their  ratios 
of  the  rescaled  data.  Then  the  regional  growth  factor  for  different  return  periods  can  be  obtained  using  equation  (11). 


(11) 


The  flood  quantiles     i  for  different  return  periods  for  the  sites  in  the  region  can  be  estimated  by  multiplying  the 
regional  growth  factors  q  (F)  by  the  scale  factor  or  site  mean     =     as  per  Equation  (11). 
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REGIONAL  GROWTH  FACTORS 

The  best  fit  distributions  for  both  the  data  are  found  as  LN3.  The  growth  factors  for  the  valley  can  be  derived  from 
the  quantile  function  of  LN3  by  applying  index  flood  procedure.  The  regional  values  of  A  ^J1,  A^and  twa  of  the  rescaled 
data  for  both  the  data  and  the  parameters  of  the  LN3  distributions  are  in  given  in  Table  7. 

The  equations  for  regional  growth  factors  derived  from  the  quantile  function  of  LN3  for  the  for  both  discharge 
and  gauge  level  are  given  in  Equation  (12),  and  Equation  (13)  respectively. 

Table  7:  Regional  L-Moment  And  Their  Ratios  Of  Index  Flood  Procedure  With  The  Parameters  Of  The  Fitted 

Distributions 


L-Moment  of  the  Rescaled 
Data 

Discharge 

Gauge 

1.0 

1.0 

Jljj    =  Trt 

0.132442 

0.017443 

0.0581456 

-0.090192 

0.152124 

0.154976 

Parameters  of  LN3  distribution 

s 

0.98606 

0.9972 

a 

0.23321 

0.03049 

K 

-0.11909 

0.18494 

DISCHARGE  DATA 

=  -0.97212  +  0.382025  (12) 
GAUGE  LEVEL  DATA 

—  =  -0.97212  +  0.382025  BDJ-«»»  *  (13) 

The  regional  growth  factors  for  different  return  periods  for  both  discharge  and  gauge  level  is  given  in  Table  8. 
The  estimated  flood  in  terms  of  discharge  and  gauge  level  for  different  return  periods  for  sites  in  the  regions  A  and  B  are 
given  in  Table  9  and  Table  10  respectively. 


Table  8:  Regional  Growth  Factors 


Data 

Growth 
Factor 

Return  Periods 

2 

5 

10 

15 

20 

30 

40 

50 

75 

100 

Discharge 

Or 

& 

0.9937 

1.0376 

1.0646 

1.0793 

1.0893 

1.1030 

1.1124 

1.1196 

1.1323 

1.1411 

Gauge 

Or 
|i 

0.9972 

1.0209 

1.0320 

1.0372 

1.0404 

1.0446 

1.0473 

1.0493 

1.0526 

1.0548 
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Table  9:  Estimation  of  Flood  Discharges  for  Different  Return  Periods 


Station 

2 

5 

10 

15 

20 

30 

40 

50 

75 

100 

1.  A.P.  Ghat 

3611.31 

4366.82 

4793.83 

5015.03 

5163.35 

5361.84 

5496.98 

5598.79 

5777.87 

5900.92 

2.Fulertal 

3908.24 

4725.87 

5188 

5427.38 

5587.9 

5802.71 

5948.96 

6059.14 

6252.95 

6386.11 

3.  Gumra 

168.26 

203.46 

223.35 

233.66 

240.57 

249.82 

256.12 

260.86 

269.2 

274.94 

4.  Dholai 

301.01 

363.98 

399.57 

418.01 

430.37 

446.92 

458.18 

466.67 

481.59 

491.85 

5.  B.P.  Ghat 

4995.28 

6040.33 

6630.99 

6936.96 

7142.12 

7416.68 

7603.6 

7744.43 

7992.14 

8162.35 

6.Fakirabazar 

395.66 

478.44 

525.22 

549.46 

565.71 

587.46 

602.26 

613.42 

633.04 

646.52 

7.  Matizuri 

1296.14 

1567.3 

1720.56 

1799.95 

1853.19 

1924.43 

1972.93 

2009.47 

2073.74 

2117.91 

Table  10:  Estimation  of  Flood  Gauge  Level  for  Different  Return  Periods 


Station 

2 

5 

10 

15 

20 

30 

40 

50 

75 

100 

1.  A.P.  Ghat 

20.54 

21.03 

21.26 

21.37 

21.43 

21.52 

21.57 

21.62 

21.68 

21.73 

2.Fulertal 

24.76 

25.35 

25.62 

25.75 

25.83 

25.94 

26 

26.05 

26.14 

26.19 

3.  Gumra 

17.6 

18.02 

18.21 

18.31 

18.36 

18.44 

18.48 

18.52 

18.58 

18.62 

4.  Dholai 

27.15 

27.8 

28.1 

28.24 

28.33 

28.44 

28.52 

28.57 

28.66 

28.72 

5.  B.P.  Ghat 

17.46 

17.88 

18.07 

18.16 

18.22 

18.29 

18.34 

18.37 

18.43 

18.47 

6.Fakirabazar 

15.83 

16.2 

16.38 

16.46 

16.51 

16.58 

16.62 

16.65 

16.7 

16.74 

7.  Matizuri 

17.17 

17.58 

17.77 

17.86 

17.92 

17.99 

18.03 

18.07 

18.13 

18.16 

8.  Matizuri 

21.75 

22.27 

22.51 

22.62 

22.69 

22.78 

22.84 

22.89 

22.96 

23.01 

ESTIMATION  OF  FLOOD  FOR  THE  UNGAUGED  SITES 

The  peak  flood  in  terms  of  discharge  and  gauge  level  at  the  ungauged  sites  of  this  valley  is  estimated  in  order  to 
give  probable  information  about  flood  which  may  be  use  in  various  engineering  works.  Multiple  regression  equation  (14)  is 
developed  by  relating  mean  discharge  flood  with  channel  width  at  the  gauging  site,  elevation  and  catchment  area  for  the 
sites  for  the  estimation  of  mean  annual  peak  discharge  flood.  The  multiple  regression  equation  (15)  is  also  developed  by 
relating  mean  gauge  level  with  channel  width  at  the  gauging  site,  and  elevation  in  order  to  estimate  mean  peak  flood  level 
for  the  sites. The  results  of  the  regression  equations  are  validated  using  chi-square  test  and  site  characteristics  of  the  one 
retained  sites.  The  critical  values  for  the  chi-square  test  between  the  observed  and  estimated  flood  quantiles  are  much  less 
than  the  critical  value  of  12.592  at  5%  significance  level  with  6  degrees  of  freedom  and  14.067  at  5%  significance  level 
with  7  degrees  of  freedom  for  both  the  data  respectively.  Hence  this  shows  that  there  is  no  significance  difference  between 
the  estimated  and  observed  values.  The  estimated  annual  peak  discharge  means  for  the  sites  using  equations  (14)  and  (15) 
are  given  in  Table  1 1  and  Table  12. 

q  =  e11 7120  (Ch)  -1  1969.  (Elev.)  3  8413.  (A) 0  8580  (14) 
A  =e°-8439.(Ch)  009l4.(Elev.)0-5620  (15) 
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Table  11:  Estimation  of  Flood  Discharges  for  Ungaged  Sites 


Station 

Channel 

(Ch) 

Elevation 

(E) 

Area(A) 

Obs.(O) 

In  (ch) 

ln(E) 

ln(A) 

ln(Q) 

Est.(E) 

2 

z  = 

£ 

l.A.P. 
Ghat 

146.4 

19.2 

18721 

3662  21 

4.9863 

2.9549 

9.8374 

8.2058 

3563  26 

2.7478 

2.Fulertal 

165.35 

31.39 

14450 

3963.33 

5.1081 

3.4465 

9.5784 

8.2848 

4073.93 

3.0026 

4.Dholai 

34.52 

33.22 

562 

305.25 

3.5415 

3.5032 

6.3315 

5.7211 

296.97 

0.2309 

5.B.P. 

Ghat 

195.24 

19.81 

24216 

5065.69 

5.2742 

2.9862 

10.0948 

8.5302 

5017.64 

0.4602 

6.Fakira 
Bazar 

46.37 

17.98 

1108 

401.24 

3.8367 

2.8893 

7.0103 

5.9946 

412.57 

0.3110 

7.  Matizuri 

55.47 

29.26 

7770 

1314.41 

4.0158 

3.3762 

8.958 

7.1811 

1325.19 

0.0876 

X2  Chi  -Square 

6.8401 

3.  Gumra 

7.56 

32.3 

2800 

170.63 

2.0229 

3.4751 

7.9374 

5.1395 

193.14 

Table  12:  Estimation  of  Flood  Gauge  Levels  for  Ungaged  Sites 


Station 

Channel 

(Ch) 

Elevation 

(E) 

Obs.(H) 

Ln(ch) 

P3 

Est 
•(H) 

2 

*  -       f  " 

l.A.P.Ghat 

19.2 

20.6 

4.9863 

2.9549 

3.0253 

19.3 

0.0870 

2.Fulertal 

31.39 

24.83 

5.1434 

3.4465 

3.2121 

25.81 

0.0375 

3.  Gumra 

32.3 

17.65 

2.0229 

3.4751 

2.8707 

19.72 

0.2179 

4.Dholai 

33.22 

27.23 

3.5415 

3.5032 

3.3043 

23.02 

0.7699 

5.B.P.Ghat 

19.81 

17.51 

5.2364 

2.9862 

2.8628 

20.1 

0.3339 

6.Karimganj 

14.02 

15.87 

4.4667 

2.6405 

2.7644 

15.43 

0.0127 

7.Fakirabazar 

17.98 

17.22 

3.8367 

2.8893 

2.8461 

16.75 

0.0132 

%2  Chi  -Square 

1.4721 

8. Matizuri 

29.26 

21.81 

4.0158 

3.3762 

3.0824 

22.38 

0.0148 
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